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Abstract. We review our evaluation of forward dispersion relations for direct fits to the different, 
and often conflicting, nn scattering experimental analyses. We find that some of the most commonly 
used data sets do not satisfy these constraints by several standard deviations. We also provide a 
consistent nn amplitude by improving a global fit to data with these dispersion relations. 

A precise knowledge of the kk scattering amplitude provides crucial tests for Chiral 
Perturbation Theory (ChPT), as well as information on light meson spectroscopy, pionic 
atom decays and CP violation in kaons. However, there are several TtK scattering data 
sets in the literature, in conflict among themselves, even within the same experiment. 
The reason is that the data are extracted from other reactions and thus with large 
theoretical and systematic uncertainties. Here we review our recent works |[lt 0, Isj] 
where we checked dispersion relations on the different sets of data and provided simple 
parameterizations of rnz scattering amplitudes consistent with such requirements. 

Fits to different sets of data 

We first consider ^^i' to data for the SO, S2, P waves, below ^1/2 ^ 1 GeV. We 
parameterize the phase shifts, 5{s), taking into account the analytic properties, zeros 
and poles of the amplitude, in a conformal expansion of cot 5 (5). 

The P wave, up to ~ 1 GeV, comes from the following fit to the pion form factor |0] 

cot5iW = ^(M2-5)|5o + 5i^^^^|; .y' = 1.05 GeV. (1) 
5o= 1-069 ±0.011, 5i =0.13±0.05, Mp= 773. 6±0.9 

{sq is the point where inelasticity begins to be nonnegligible). The fit is seen in Fig. la. 

For the S2 wave at low energies, we first fix the Adler zero at zi = and fit only the 
low energy data, s^^^ < 1.0 GeV; later on we allow Z2 to vary. We have 

cotSfW = ^^|so + S,^^;^}. .„"^= 1.05 GeV. (2) 

2k s-2z2 { VS + \/so-s J 

5o= -80.4±2.8, 5i = -73.6±12.6; (3) 
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FIGURE 1. a) P wave TTTT phase shifts versus Eq.Q (solid line below 1 GeV), which is not a fit to 

these data, but to the pion form factor |01- The uncertainty is the hne thickness. Above 1 GeV, the dotted 
line and error are as in f3\. b) SO phase shifts and error band as given by Eq.0 below 1 GeV, and from 
j3| above. The K14 and K2:^ decay data are not shown, (see our L3J for details), c) Continuous line: The fit 
to / = 2, D-wave phase shift data. Broken line: fit improved with dispersion relations. The experimental 
points are from |8]. 
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Except for the very reliable K14 and K ^ 2n decay experiments |7], that we always 
fit, the SO data is very confusing and thus we have adopted two approaches. In the first 
method, the ''global fit", we fit averaged phase shift data between 0.81 GeV < 5^/^ < 
0.97 GeV, (where the experiments agree within 1.5 o), composing their errors carefully. 
The fit, shown in Fig.l.b, is valid for ^^/^ < 0.95 GeV, and with the Adler zero fixed at 
zo — Mji, corresponds to 
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3; = 2 1.04 ±0.70, x = 0±2.6. 



The second method is to fit only Ki^ and K data, or to add to them the data from 
the various experimental analyses separately. The results can be found in Table 1. 

The D2 inelasticity is negligible below 4A/p. A pole term is needed since the data 
gives a small negative phase above ~ 500 MeV, but the Froissart-Gribov representa- 
tion, yields |9] a positive scattering length, a^2^ = (2.72 ±0.36) x lO^^M^^, which is 
included in the fit. Also, the inflection seen in data around 1 GeV asks for a third order 
conformal expansion. So we write 
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(23.7 ±3.8) X 10^, 



And we find Bq = (2.4 ± 0.3) x 10^ Bi = (7.8 ± 0.8) x 10^ B2 
A = 196 ± 20 MeV. The fit is shown in Fig. 1 .c. 

For brevity, we do not discuss here the DO and F waves as they do not present special 
features, nor the intermediate energy region 1 GeV > > 1.42 GeV, but a detailed 
account of their parameterizations can be found in [3]. 

We also need the imaginary part of the scattering amplitude at 5^/^ > 1.42 GeV, that 
we take from a Regge fit to experimental data ( see [2] for details and the slightly 
improved rho residue of ||3|]), shown in Fig.2. As discussed in 10,0] > standard Regge 
factorization describes experiment [113,0] and is consistent with crossing sum rules. 
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FIGURE 2 The 7r;r cross sections. Experimental points from flQ^ . The stars at 1.38 and 1.42 GeV (PY) 
are from the phase shift analysis of experimental data given in 0|. Continuous lines, from 1.42 GeV 
(PY): Regge formula, with parameters as in fl] (the gray bands cover the uncertainties). Below 2 GeV, 
the dotted line corresponds to the n^n^ cross section from the Cern-Munich analysis of Hyams et al.|^ 
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TABLE 1 Fits to K decays and different data sets of SO phase shifts. PY is our global fit, Eq.0, and its 
Bq and Z?i uncorrelated uncertainties can be obtained from Eq.@. Note that many of these sets have very 
large x^/d.o. f. for the It = 1 and n'^n'^ dispersion relations below 950MeV. Let us also remark that some 
of them (like that of Sol. E) have a relatively small x^/d-o. f. just because it has huge uncertainties in its 
parameters. Ideally one would require a low x^/d.o.f. with small uncertainties in the parameters. 



Checking and improving amplitudes with forward dispersion relations 

Let us study how well the previous fits to data satisfy three independent scattering 
amplitudes. We choose the ?-symmetric or antisymmetric combinations, that form a 
complete set: Fqq = F{'nP'nP Tt^n^), Fo+ = F{7t^7t~^ Tt^Tt^), and the t channel 
isospin one amplitude, F^'''"^\ Final uncertainties are small for the two first, since they 
depend only on two isospin states, and their imaginary parts are sums of positive terms. 
Thus, we find two dispersion relations by choosing either F = Fqq or F = Fo+ in 

. s{s-4MI) r / (2s' -4-Ml)lmF(s') 

ReF{s)-F{4Ml) = ^ ^P.P. / d/ ,^ , \, , ^ , ^ ' ^, (5) 



By setting s = 2M^, and F = Fqq, we find a sum rule important to fix the Adler zeros. 



Finally, for isospin unit exchange, which does not require subtractions, 

ReF(^'-i) 5,0 = ^ P.P. / ds'—, — ^ (7) 

at threshold this is known as the Olsson sum rule. 

Depending on the method we use to fit the SO wave we find the results in Table 1, 
where, we have separated on top those fits to data with a total /d.o.f. < 6 for the 
n^n^ and = 1 dispersion relations up to 0.925 GeV, a fairly reasonable /d.o.f. 
since these fits were obtained independently of the dispersive approach. 

However, in Table 1 we also list the very frequently used t and ^-channel solutions 
of Estabrooks and Martin those of Protopopescu et a/.j^], from Table VI, VIII and 
table XII, as well as the solution A of Grayer et al. 16] . Their It = I plus n^n^ dispersion 
relation total x^/d-o.f. is surprisingly poor: 11.3, 10.1, 7, 6, 7.5, 9.9, respectively. 
Therefore, any result that relies heavily on these sets should be taken very cautiously. 

We have also improved the previous low energy fits parameters by fitting also the 
dispersion relations up to 0.925 GeV, thus obtaining parameterizations more compatible 
with analyticity and s — u crossing. Improving from Eq.©, we find, in M^^ units. 



S2; 5^/2 < 1.0 
P; 5^/2 < 1 05 
D2; 5^/2 < 1.42 



SO; 5^/2 < IniK : Bq = 17.4±0.5; Bi = 4.3 ± 1.4; 

= 790 ±21 MeV; zo = 195 MeV [Fixed]; 
Bo = -80.8 ± 1.7; Bi = -77 ± 5; Z2 = 147 MeV [Fixed]; 
Bo = 1 .064 ± 0. 1 1 ; 5i = 0. 170 ± 0.040; Mp = 773.6 ± 0.9 MeV; 
Bo = (2.9±0.2) X 10^; Bi = (7.3±0.8) x 10^; 
52 = (25.4±3.6) X 10^; A = 212±19. (8) 

The DO and F waves do not change appreciably and we just refer to [0] for details. In 
Fig. 3 we show the improved curves for SO and S2, and that of D2 in Fig. 1 .c. 

Concerning the improved fits to individual data sets, we get somewhat different results 
for SO, listed in Table 2 with the x^/d.o.f. of each forward dispersion relation and the 
standard deviations for the sum rule in Eq.® (which are more than four for K decay plus 
the Grayer B or E or Kaminski improved solutions). For other waves, no matter what SO 
fit is used, we find very similar values to those in Eq.®. This can be checked in Fig.3.b, 
where we show the improved "K decay + Grayer Sol. B" S2 wave. Even though it is the 
one for which we obtained the most different central values for the SO wave compared 
with those given in Eq.® , it falls within the uncertainty of our improved solution. 

In summary, with forward dispersion relations, we have checked the consistency of 
different phase shift analyses available in the literature. Surprisingly, some of the most 
frequently used phase shift sets do not satisfy these dispersive constraints and sum rules, 
and should therefore be used cautiously. We have provided a simple parameterization 
of nn scattering consistent simultaneously with some data sets and all three forward 
dispersion relations, that we hope could be of use for future studies of scattering. 
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FIGURE 3. a) The improved SO phase shift (PY improved, EqEJ, the global fit (PY from data, the SO 
in Eq. Q), and the improved solutions "K decay only" and "Grayer C" of Table 2 (almost on top of PY 
improved), b) S2 improved Phase shift (PY improved, Eq. (|8}); global fit (PY from data, Eq. 0) and the 
improved parameterization with K decays and So. B of Grayer et al.|6]. 



TABLE 2. Improved fits (only those already with /dof < 6 in Table 1). Although errors are given 
for the Adler zero, we fix it when evaluating other errors, to break the otherwise very large correlations. 
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